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Abstract 

Zero temperature states of matter are holographically described by a spacetime with 
an asymptotic electric flux. This flux can be sourced either by explicit charged matter 
fields in the bulk, by an extremal black hole horizon, or by a combination of the two. 
We refer to these as mesonic, fully fractionalized and partially fractionalized phases 
of matter, respectively. By coupling a charged fluid of fermions to an asymptotically 
AdS^ Einstein-Maxwell-dilaton theory, we exhibit quantum phase transitions between 
all three of these phases. The onset of fractionalization can be either a first order or 
continuous phase transition. In the latter case, at the quantum critical point the theory 
displays an emergent Lifshitz scaling symmetry in the IR. 



1 Context 



A field theory with a global U(l) symmetry may be placed at a finite charge density (</'). 
It is of interest to characterize the resulting possible low energy, zero temperature phases of 
matter. In this paper we will consider phases of matter where no symmetries (i.e. neither the 
global U(l) nor the spacetime symmetries) are spontaneously broken. Lorentz invariance 
will of course be broken by the charge density itself. 

The holographic correspondence [1] gives a set of theories for which this question may be 
posed nonperturbatively. The charge density is dually described by an asymptotic electric 
flux *F in the holographically dual spacetime. This electric flux must be sourced in the 
interior of the spacetime. In the classical gravity limit there are two qualitatively distinct 
sources for this flux: charged matter in the bulk or a charged event horizon (see e.g. [2]). 
The sharpness of this distinction is made possible by the fact that the large N limit of 
theories with gravity duals necessarily introduces a parametric hierarchy between a few 
'light' bulk fields and a large number of 'black hole microstates', see e.g. [3J. 

Charged matter fields in the bulk correspond to (composite) gauge invariant charged 
operators in the dual field theory. At a classical level in the bulk, such fields must be 
fermionic if they are to contribute to the charge density without condensing or breaking 
a spacetime symmetry. It has recently been understood that the charge density carried 
by such 'mesinos' (i.e. fermionic meson-like operators) retains key features of conventional 
Fermi liquids. In particular, the sum over the volumes of all Fermi surfaces in the mesino 
correlators equals the total charge density carried by the mesinos, in accordance with the 
Luttinger theorem [H [5j |6] . The phenomenology of these Fermi surfaces can be unconven- 
tional due to dissipation into a locally quantum critical sector El EJ [TUl Ell ECS H3j or 
due to the existence of many very closely spaced Fermi surfaces [21 HI [15]. Nonetheless, 
all of these features can be understood without holography as conventional Fermi surfaces 
coupled to a critical sector in a way that allows the fermion self energy diagrams to be 
reliably resummed at large N \16\ [T7] . 

Charge emanating from behind an event horizon, in contrast, appears as a mismatch in 
the total Luttinger count of gauge invariant Fermi surfaces of 'light' operators [2]. It has 
been argued that such a mismatch is a defining characteristic of fractionalized phases of 
matter, in which the remaining charge density is thought of as being carried by the gauge- 
charged 'quarks' [18J. Indeed, by analogy with the holographic association of deconfined 
phases with horizons |19j it is natural to associate flux emanating from a horizon with 
fractionalized charge carriers [H [2] . The association of charged horizons to fractionaliza- 



1 



tion remains to be made as precise as that of horizons in general to deconfinement, as we 
currently lack an order parameter for charge fractionalization analogous to the Polyakov 
loop. 

The onset of fractionalization may be an important ingredient in understanding contin- 
uous quantum phase transitions in which a Fermi surface undergoes a major re-organization 
or disappears entirely \20\ \21\ [22] . Such transitions may be important in turn for elucidating 
the puzzling phenomenology of, among other materials, the heavy fermion compounds and 
cuprates. However, conventional field theoretic treatments of gauged Fermi surfaces face 
serious technical challenges, especially in 2+1 dimensions [23 1 . 

In this paper we present and explore a zero temperature holographic framework in 
which a fraction of the bulk electric flux emanates from behind a horizon, while a fraction is 
sourced by mesinos. The extent of the fractionalization - i.e. flux from the horizon - can be 
dialed using a relevant operator in the dual UV field theory. For large positive couplings, 
all the charge is fractionalized. At large negative couplings, all the charge is accounted 
for by mesino Fermi surfaces. For a finite range of intermediate couplings, both a charged 
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Figure 1: Phase diagram. By tuning a relevant coupling in the relativistic UV theory at 
fixed nonzero chemical potential, the charge density ranges from fully fractionalized to fully 
mesonic. The fractionalization transition is controlled by an IR fixed point with dynamical 
critical exponent z > 1. The IR critical point is sometimes preempted by a first order phase 
transition. 
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horizon and mesinos are present. This range of coupling is bounded on the positive end by a 
continuous phase transition and on the negative end by a continuous or first order transition, 
depending on the parameters of the theory. The continuous transition is controlled by a 
relevant operator in an emergent IR Lifshitz fixed point. This phase diagram, illustrated 
in figure [T] above, is the main result of this paper. The three phases that appear have 
important similarities with the FL, FL* and NFL phases discussed in |18j . 

2 Bulk theory and equations of motion 

To discuss fractionalization holographically, the minimal bulk ingredients we need are the 
metric, a Maxwell field and some charged fermionic matter, dual to the mesino operators. 
The backreaction of the bulk density of fermions on the metric and electromagnetic fields 
is important [24]. For this purpose, it is technically convenient, and perhaps natural, to 
consider the fermions in a limit in which they admit a coarse-grained description as an ideal 
fluid. In this WKB limit, the fermion wavefunctions become highly localised - more detailed 
discussions of the regime of validity of the fluid description can be found in |25l [2"6"l |2"71 H] . 
It is possible to go beyond the fluid limit using the techniques recently considered in [6], 
but we shall not do so here. 

Let us briefly review the minimal Einstein-Maxwell-charged fluid framework. There are 
two types of finite density, zero temperature solutions. If the fermion mass is sufficiently 
large relative to its charge, then the chemical potential in the bulk is insufficient to populate 
a Fermi sea and the bulk solution is an extremal Reissner-Nordstrom black hole. Following 
the comments in the introduction, we think of this as a fully fractionalized phase. These 
extremal black holes are stable against pair production of fermions outside the horizon - 
because the fermions have more mass than charge (in an appropriate dimensionless sense), 
pair produced fermions fall through the event horizon rather than being pushed in opposite 
directions and discharging the black hole. If the fermion is sufficiently light, however, then 
a Fermi sea is populated in the bulk and, solving the equations of motion, it is found that 
all the charge resides in the fermion fluid, with no charge emanating from a horizon [26] . 
In this case we have a fully mesonic phase that has been dubbed an 'electron star'. 

Upon heating up the electron star, charge is gradually transferred from the fermion 
fluid to the horizon. At a critical temperature there is a third order phase transition, above 
which all the electric flux is sourced by the horizon [281 [29] . This can be thought of as a 
very mild cousin of the Oppenheimer-Volkov bound on the mass of neutron stars, which 
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more typically give first order phase transitions |25[ [27] . This is the kind of fractionalization 
transition we will be interested in. However, we would like the transition to occur at zero 
temperature, where Fermi surfaces are more sharply defined, as a function of some coupling 
in the field theory. For this we need more ingredients. 

A commonly considered generalization of Einstein-Maxwell theory is to add a dilaton 
field. We will consider the following Einstein-Maxwell-dilaton-charged fluid action 

z = 2^2 ( R ~ 2 ( v$ ) 2 " ^W) ~ Z ^r F ^ Fab + p0*O • (i) 

Here p is the pressure of the fluid. For our irrotational, zero temperature fluids we can take 
the local chemical potential in this action to be 

tf oc .=g ab A a A b , (2) 

see for instance [26], following [30] . 

Without the charged fluid, this type of action has been considered recently as a model 
for QCD at finite baryon density |31j . as well as to produce a wide range of extremal near 
horizon geometries with a view to phenomenological modeling of frequency and temperature 
scaling laws in condensed matter systems (32J, [331 EI] • Various interesting observations about 
fermionic probes of Einstein-Maxwell-dilaton backgrounds were made in |35[ 1361 [37] , but 
will not be directly related to our discussion. Within the context of QCD, the dilaton is 
a simple way to incorporate a relevant operator that can mimic the effects of the running 
QCD coupling and the ultimate onset of confinement, see e.g. [381 E2]- We shall use the 
dilaton in a similar way, choosing its mass so that it is dual to a relevant operator. 

We will study solutions to this theory with metric and Maxwell potential 

ds 2 = L 2 (-f{r)dt 2 + g(r)dr 2 + & , (3) 
A = —h(r)dt. (4) 
Thus in particular the local chemical potential 

Woe. = Tf- (5 

The UV boundary will be at r = 0, while the IR interior will be at r — > oo. The energy 
density, charge density and pressure of the fluid are not independent variables, but are 
rather determined at each point by this local chemical potential. If we introduce the rescaled 
dimensionless variables 



then the equation of state for the zero temperature Dirac fermion of mass m is 

h h 

p = /3 e 2 \/e 2 — m 2 de , a = ft I e\f e 2 — rh 2 de , —p = P Ti^ ■> CO 

Jrh Jrh V J 



where 

Jr2 1 
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Note that these local thermodynamic quantities obey (as they must) 

#' = (£)'»■ w 

A more leisurely exposition of the above may be found in [26]. 

On the ansatz just described, the equations of motion following from our action are 
found to be 

1 (f , 9' , A , gh , , - 



, , +- + - + ^=«t + 2$" = 0, (10) 



2\ f g rj 4 4/ 



d 



\r 2 Vfg) r 



a = 0. (13) 



dr 

These equations are seen to imply the following additional conservation equation, cf. 

d ( 2r 2 Z(<S>)hh>-{r 2 f)> \_ 

dr\ r*Vf9 J ' { ' 

Such additional conservation equations are important for the first law of thermodynamics 
to be satisfied, as they connect data at the horizon and at the asymptotic boundary. For the 
zero temperature spacetimes to be considered in this paper, this charge will in fact vanish, 
so that 

2r 2 Z{<$>)hti - (r 2 f)' = 0. (15) 
This can be checked explicitly on the all the near horizon geometries we present below. 

3 Critical points and near horizon solutions 

3.1 Asymptotic AdS^ solution and thermodynamics 

We will choose the dilaton potential such that there is an AdS± solution to our equations of 
motion, with no fermion fluid. We will use this solution to set the UV boundary conditions 
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on our theory. The general near boundary, r — > 0, expansion will be characterized by several 
constants. These are the boundary speed of light c, the energy density E, the source and 
expectation values 4>q and (O) of the operator O dual to the bulk dilaton $, and finally the 
charge density Q and boundary chemical potential £i. The various hats indicate that we are 
dealing with quantities that have been rescaled by factors of e, L, k. At this point we will 
assume, for ease of computation, that the potential has the form near <3? = of 

V($) = -6 - 4$ 2 + 0($ 4 ) , (16) 

so that the dilaton has mass squared of —2/L 2 . This is above the AdS± Breitenlohner- 
Freedman bound; hence the dual operator O has scaling dimension A = 2 in the UV 
theory. We will also set Z(0) = 1 without loss of generality. The near boundary expansion 
then takes the form 

/ = ^(l-(^+^o(0))r 3 + ---) , (17) 
g = I(l-^ r 2+^-0 o (C»))r 3 + ---) , (18) 
h = c (# - &■ + ■■■) , (19) 
$ = 0or + ^r 2 + ■■■ . (20) 

Thus the theory is asymptotically characterized by three sources {c, (po,fi} and three ex- 
pectation values {E, (O), Q}. 

Using the equations of motion, the bulk Lagrangian ([T]) is found to become a total 
derivative on shell 

r— c \ L 2 d ( 2Z^)hh'-f \_L 2 d f 

V ^Ion-shell K 2 dr y 2r 2 ^/Jg~ J k 2 dv v^yj~fg K ' 



In the second equality we used the conservation law (14). To obtain the free energy from 
the on shell action, we need to add the usual boundary counterterms, see e.g. [ID] and 
references therein, 

£ «=M*-t)-i* 2 - (22) 

Here K = j^V^riu with 7 the induced metric on the asymptotic boundary and n an 
outward pointing unit normal. The free energy density is now found to be 

il 1 lim ( / ^Cdr + V^£ c .t.L ] (23) 



V C e->0 

= E-fiQ. (24) 
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Here we have used the above expansions and imposed that we are at zero temperature and 
hence there is no contribution from the endpoint of the radial integral in the interior of 



the spacetime. Evaluating the conservation law (14) in the asymptotic region furthermore 
implies the relationship 

^E-fiQ + ^MO) = 0. (25) 

This expression follows from tracelessness of the energy momentum tensor for a conformally 
invariant theory, E = 2P, together with the effect of the relevant deformation <j)Q. Consider 
changing the volume at fixed chemical potential and zero temperature. Differentiate the 



free energy (24) and simplify using the first law of thermodynamics, d(EV) = —PdV — 
(0)d(<f)oV) + fid(QV). Extensivity of the free energy allows us to integrate back again to 
give £l/V = -P - {0)d(4>oV)/dV = -P - \<j) Q {0). In the second step we used the fact 



that 4>q has scaling dimension one. Using the tracefree condition, we obtain (25) 



3.2 Lifshitz fixed point 

In the presence of the charged fluid, our system admits a second scale invariant solution. 
This will be a Lifshitz spacetime |41j characterized by a dynamical critical exponent z. The 
dilaton will be constant, and therefore the physics governing these solutions is essentially 
the same as that of the IR geometry of electron stars without the additional dilaton field 
\24\ [26] . The difference will be the presence of a relevant perturbation of the Lifshitz fixed 
point. The system must therefore be tuned in the UV to hit this IR Lifshitz solution. 
We will discuss the renormalization group trajectories below, as well as the nature of the 
continuous phase transition that occurs when the system is driven through the Lifshitz fixed 
point. 

The Lifshitz solution takes the form 

, 9 r 9 / dt 2 dr 2 dx 2 + dy 2 \ , eL dt _ . , 

ds = L \-^-z + 9l^t + 5-M , A = h L -, $ = </>L- 26 

So there are four constants to solve for: {z,gL,hi,,(f>L}- The dilaton equation of motion 



(12) implies that the field must be at the minimum of an effective potential receiving 



contributions from the potential and also from the electric flux 

Kff.(^) = QlV'^l) - z 2 h 2 L Z'(4> L ) = . (27) 



In general this equation may have multiple solutions. Meanwhile, the Gauss law (13) 
selfconsistently relates the potential hi to the constant charge density ct(^) 

2zh L Z{cp L ) = g L a{h L ) . (28) 
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The conservation law ( 14 ) gives a simple expression for the dynamical critical exponent 

1 



(29) 



1 - h\Z{$ L ) ' 

Thus in particular z > 1 if Z(<&) is positive. The remaining equation of motion is then seen 
to give, after simplification using the above expressions, 

(1 + z)(2 + z) = 2g L p(h L ) - g L V(^ L ) . (30) 

For a given theory with V(&) and Z($), we can use these equations to eliminate z and gz 
and obtain coupled equations for {hL,(f>i} that must be solved numerically. 

For our Lifshitz spacetimes, the conserved electric flux J K2 * [Z($)F] ~ r -2 — > in the 
interior as r — > oo. Thus there is no flux emanating from the Lifshitz 'horizon', and the 
Lifshitz IR describes a fully mesonic configuration. 

At this point it is convenient to fix our model. For the remainder of the paper, we 
choose the following simple form for the scalar potential and coupling to the Maxwell term 

Z($) = e 2$/ ^, F($) = -6cosh(2$/ v / 3) • (31) 

This form was largely chosen for convenience of numerical implementation. The potential 



satisfies our assumed expansion (16) at small fields, so that the operator O dual to the 
dilaton is relevant in the UV with dimension A = 2. The fact that Z is not symmetric 
in $ implies that = will not be a solution to the equations of motion once there is 
electric flux, even when the dual operator O is not sourced explicitly. Many features we 
are interested in will not depend on the details of Z and V, we shall comment on this as 
appropriate belowQ 



With the expression (31) it is straightforward to solve the algebraic equations for 
{z, gz, <Pl} for a given {/3,m}, characterizing the fermion fluid. The Lifshitz solution 
can only exist if the local chemical potential is large enough to populate the Fermi sea, 
which requires 

(32) 

The above equality may be numerically solved to find m max . • The following figure [2] shows 
m max , as a function of (3. 



1 A positive energy theorem for the Einstein-dilaton sector of our theory is equivalent to the existence of a 
superpotential W [42], which with our normalization must satisfy V = 2(d<s>W) 2 — 6W 2 . While we have not 



found W analytically for our potential (31 1, we have checked numerically that it exists and has asymptotic, 
large <&, behavior W ~ cosh 
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Figure 2: Maximal mass m max . as a function of j3. The Lifshitz solution requires rh < rh a 



Given the Lifshitz solution, we need to characterize the deformations away from the 
fixed point. To this end, write 



h 



hL 
1 



1 + 5hr 



[1 + Sfr 



9 - 



9L 



(1 + 



<r M ) , 



(33) 
(34) 



Because the total order of our equations of motion ftWn - (13) is six, we can expect to 



obtain six different possible exponents M. These exponents will consist of three pairs of 
sources with momentum scaling dimensions [go] = M_ and expectation values with scaling 
dimensions [(C)] = M + satisfying 



M + + M_ 



2 + z. 



(35) 



This relation comes from the fact that each pair of modes describes the insertion of an 
operator f dtd 2 x go O at the Lifshitz fixed point. Here O represents any one of the three 
operators in the Lifshitz IR theory, not the operator dual to the dilaton. 

A key fact we will wish to characterize is whether the operators that can perturb the 
Lifshitz fixed point are relevant or irrelevant. There is a universal deformation with [go] = 
and associated expectation value with dimension 2 + z that corresponds to placing the 
theory at a finite temperature (cf. for instance [431 [33], I26j). We will not want to source 
this deformation. We then find that of the remaining two operators, one is irrelevant and 
one is relevant. The irrelevant operator has [go] < 0, and hence becomes more important 
towards the UV boundary at r — > 0, while the relevant operator has [go] > 0. 



9 



Consider the irrelevant deformation first. We will find below, when we consider full RG 
flows, that the irrelevant operator allows us to flow up from the Lifshitz fixed point to the 
relativistic UV fixed point. See figure [I] above. In figure [3] we show the dimension of the 
irrelevant operator, and also the value of the dynamical critical exponent, for a range of 
allowed values of {f3,m}. The dynamical critical exponent and dimension of the irrelevant 
0.8 F 




Figure 3: Contour plot of the dynamical critical exponent z (left) and dimension [go] of the 
irrelevant operator (right) as a function of {/3,m}. 



operator will determine the scaling of quantities as the Lifshitz point is approached from the 
UV. In particular, they will determine the low temperature scalings of observables. It is of 
interest to heat up the phase diagram in figure [TJ to describe the thermodynamic crossovers 
around the critical point, along the lines of the recent interesting work |44j . 

The presence of the relevant deformation means that we will need to tune the UV theory 
if we want to hit the Lifshitz fixed point in the IR. Indeed, we have a dimensionless ratio 



of relevant couplings that can be tuned in the UV: 0o/A from section 3.1 above. We can 
expect that a critical value of this ratio will correspond to the critical red line in figure [I] 
that flows between the two fixed points. A little away from the critical ratio, as the flow 
runs close to the IR Lifshitz fixed point, it will pick up the relevant coupling in the IR and 
flow away from the Lifshitz fixed point. We can therefore anticipate that the Lifshitz fixed 
point mediates a transition between two distinct phases of matter. Indeed we will find below 
that these correspond to mesonic and fractionalized phases. One important subtlety here is 
that the dimension of the relevant operator turns out to be complex for a range of values of 
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{/§, m}, as we show in figure [4] below. This phenomenon has been observed previously about 
Lifshitz fixed points, e.g. in [35]. The two powers M± are complex conjugates of each other. 
In this case (35) necessarily implies Re [go] = (2 + z)/2. Complex scaling exponents are 




Figure 4: Contour plot of the real (left) and imaginary (right) parts of the scaling dimension 
[go] °f the relevant operator as a function of {/3,m}. A nonzero imaginary part indicates 
an instability of the Lifshitz solution, leading to a first order phase transistion. 



known to indicate the presence of an instability, see [H] and [36j for AdS-z x K 2 and Lifshitz 
examples, respectively. Indeed, we will find below that when the scaling dimensions are 
complex, the Lifshitz solution is unstable and the continuous phase transition is preempted 
by a first order transition between mesonic and partially fractionalized phases. 

At the Lifshitz fixed point we have just described, the dilaton was stabilized in the IR 
at the value $ = In the following two subsections we will see that the two phases on 
either side of the Lifshitz fixed point are characterized by runaway IR behaviors in opposite 
directions: $ —> ±00. This distinction is possible because we chose Z in (31) not to be 
symmetric under $ -H- — 



3.3 Near horizon geometries with flux (fractionalized) 

In this subsection we look for IR geometries, i.e. r — > 00, such that the dilaton becomes 
large and positive. It follows from (31) that Z — > 00, and hence the effective Maxwell 



coupling is becoming small. The solution can be determined as a series expansion together 
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with a perturbation that will allow the flow to be integrated up to the UV: 




71=1 



OO 



n=l 

oo 



oo 

$ = ^\ogr + Y,% + ^r N ■ (39) 

n=l 

The coefficients {f n ,9n, h n ,p n } are uniquely determined by the equations of motion, while 
the perturbation {8f, 5g, Sh, 6(f)} has an overall free magnitude. This irrelevant pertubation 
has 

N = 2- 2^19/3 < 0. (40) 

We can expect that different magnitudes for the perturbation will map onto different values 
of the dimensionless ratio of relevant couplings </>o/A in the UV theory. 

From the IR expansion above we have that at leading order the local chemical potential 
Aloe. ~ - > 0. It follows that for any nonzero fermion mass m, there will be no fluid in 
the far IR of the spacetime. Because (i\ oc . is increasing as we go towards the UV, it may be 
the case that there is a fluid of fermions in an intermediate region of the spacetime, similarly 
to the finite temperature spacetimes considered in |28[ [29] . We will see this explicitly later 
when we perform a full integration of the equations. 

The conserved electric flux, in contrast, tends to a constant as r — > oo, leading to a 
nonzero flux emanating from the null singularity in the far IR: lim L 2 * [Z(&)F] ~ const.. 
The emergence of flux from behind a 'horizon' indicates that these IR geometries describe 
phases of matter that are at least partially fractionalized. 



The IR solution (36) - (39) is not scale invariant as r — > oo. This is manifested in the 
fact that g dr 2 scales, as well as the logarithmic dependence of the dilaton. The curvature 
invariants are not constant and indeed become large in the IR. Lack of scale invariance 



implies that we have a set a dimensionful scale to unity in (36) - (39). We shall only 
compute dimensionless ratios of physical quantities in the UV, and therefore this scale will 
drop out. While the effective 'string coupling' is becoming small in the limit Z — > oo, the 
divergence of curvature invariants indicates that higher derivative corrections eventually 
become important at the lowest energy scales. 



12 



The IR geometry we have just presented, without a charged fluid, is a special case of 
the families of IR solutions discussed in, e.g. [Ml ESJ [36] • Those papers also considered 
the effects of a small temperature on the geometry, by finding near-extremal black hole 
solutions, and computed the temperature dependence of the specific heat. Specialized to 
our case, their formulae show that c s ~ T at low temperatures. A power law dependence 
of the specific heat indicates that this phase of the system contains gapless excitations. 

3.4 Near horizon geometries without flux (mesonic) 

In this subsection we look for IR geometries such that the dilaton becomes large and neg- 



ative. It now follows from (31) that Z — > 0, and hence the effective Maxwell coupling is 
becoming large. We found the following IR series expansion and irrelevant perturbation in 
this case: 

/ = ^ll + r,d^ + 6fr P ) , (41) 

+ Sgr p j , (42) 
5hr p - 2 ^ , (43) 

n=l 

Here fg^o = ft{ho), there is a nonvanishing charge density of fluid in this case, and the 
{fn, 9n, h n ,p n } are similarly uniquely determined by the equations of motion. As previously, 
the perturbation {5f, 5g, 6h, 5(f)} has a free overall magnitude, and generates an RG flow to 
the relativistic UV fixed point. The exponent 

p ='-!)/ i+ 4(S)< - (45 > 

In contrast to the previous subsection, here we see that fi\ oc , — > Hq > rh, so that a 
fermion fluid is present in the far IR. Also in contrast to the previous section, here we find 

lim r„ 2 * [Z($>)F] ~ r~ 7 / 3 — > 0, so that there is no conserved flux emanating from the 

r— too 

'horizon' in the far IR. These solutions therefore describe a fully mesonic phase, with all 




flux sourced by the fermion fluid. In fact, the solution (41) - (44) can be understood in 
the following simple way. The leading behavior for the metric and dilaton, i.e. {/, g, $}, is 
precisely that of the 'domain wall' spacetime with no electric flux or fluid, i.e. h = a = 0. 
We proceed to sprinkle a density of charged fluid onto this domain wall spacetime. These 
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fermions then self-consistently source an electric field according to the Gauss law in the 
domain wall background, giving at leading order in the IR the expression that we have 
already quoted 

^Uo = Hho) ■ (46) 

Because Z — > in the far IR, we might expect that the backreaction of the electric flux 
and fermion fluid on the spacetime and dilaton will be subleading as r — > oo. Indeed, this 
is what occurs. 

As in the previous subsection, the far IR solution is not scale invariant and suffers from 
divergent curvature invariants leading to a null singularity. In the present case, the 'string 
coupling' is also divergent in the IR, so string loop effects will become important at the 
lowest energy scales. 

Turning on a very small temperature, it is plausible that the specific heat will be that of 
the finite temperature domain wall solution, with the effects of the electric flux subleading 
in the IR, as we found for the zero temperature solution. That problem was studied in some 



generality in, for instance, [33]. For our choice of potential V in (31), their results show 
that any nonzero temperature leads to a black hole solution, with specific heat c s ~ T 3 / 2 at 
small temperatures. The system is therefore again gapless, although there are fewer degrees 
of freedom at the lowest energy scales than we found in the fractionalized phase. 



By a slightly different choice of potential V in (31) we can easily find a model in which 
the domain wall spacetime, and hence presumably the mesonic phase we are presently 
discussing, is gapped. In this case a black hole does not form until the system is heated 
above a deconfinement temperature Tc, see e.g. [33]. While it is natural to associate 
fractionalized phases with deconfinement, the two are logically distinct. In the model we 
have chosen to focus on, the 'glue' remains deconfined throughout, and yet the system 
exhibits a fractionalization transition in the charged sector. The situation is similar to 
transitions in some holographic models of QCD, wherein the quarks remain bound into 
mesons above the deconfinement temperature, and only 'melt' at a higher temperature [47] . 



4 The full solutions and phase diagram 

Putting together all the ingredients above we can map out the zero temperature phase 
diagram of our model. In the UV theory we have two relevant couplings, the chemical 
potential ft and the coupling 0o for the operator O dual to the dilaton in the bulk. We 
can therefore characterize the theory as a function of the dimensionless ratio (f>o/jl. In 
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the following we present results for the free energy f2, charge density Q and the ratio of 
fractionalized to total charge Qtac./Q as a function of 4>o/£i. 



The equations ( 10 ) - ( 13 ) are solved by starting in the far IR with the metrics of sections 
3.3 or 3.4 The only free parameter is the strength of the irrelevant perturbation about the 



IR solutions. For a given value of this parameter we numerically integrate the solutions up 
to the asymptotic region at small values of r. Near the boundary we can read off the ratio 
</>o /A as wen as thermodynamic quantities, by fitting to the boundary expansions of section 



3.1. In performing the integration we must set the fluid variables {a,p,p} to zero whenever 



Aloe. ( r ) < fa and keep them in the equations otherwise. 

In figure [5] below we plot the free energy as a function of the relevant coupling 4>o at 
fixed chemical potential p.. We always plot dimensionless ratios. Two cases are plotted, 
corresponding to different values of the fermion fluid parameters {m,/3}. The properties 
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Figure 5: Free energy as a function of the relevant coupling 4>q for {m,/3} = {0.1,20} 
(left) and {fa, $} = {0.5, 10} (right). Dashed blue lines indicate a mesonic phase while the 
solid red lines indicate a partially fractionalized phase. The black dot denotes the location 
of the Lifshitz fixed point. The left plot exhibits a continuous fractionalization transition, 
while the right plot exhibits a first order transition. 



of the Lifshitz fixed point for the two values we chose are shown in the following table. 
These values are illustrative of continuous and first order fractionalization transitions and 



{fa J} 


z 


[ffirrel.] 


[diet] 


stability 


{0.5,10} 


2.62 


-2.43 


2.31 ± 1.25 1 


unstable 


{0.1,20} 


1.42 


-1.11 


1.52 


stable 



Table 1: An example of a stable and unstable Lifshitz fixed point. 
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can also be read off from the plots in section 3.2 above. As we see in the free energy plots 
of figure [5j when the relevant operator about the Lifshitz fixed point has a complex scaling 
dimension, then the Lifshitz point is in fact a local maximum rather than a minimum 
of the free energy. The continuous fractionalization transition in this case occurs on an 
unstable branch of solutions, and so the continuous transition is preempted by a first order 
transition, as the two stable branches (the minima) cross. In these plots and throughout 
this section, we use dashed blue lines to denote mesonic phases and solid red lines to denote 
phases with partial or total fractionalization. Points on the blue lines are obtained by 



integrating outwards from the near horizon geometry of section 3.4 while points on the red 



lines are obtained by integrating out from the near horizon geometry of section 3.3 As an 
independent check of the numerics, we also integrate out from the Lifshitz solution. The 
plots in figure [5] are very zoomed in to the critical point; it is necessary to use high precision 
numerical integration and to read off the asymptotic quantities with some care. 

Combining the information in figures [3] and [4] we can conclude that continuous fraction- 
alization transitions, with real scaling dimension for the relevant operator at the Lifshitz 
fixed point, can only occur for either relatively small or relatively large dynamical critical 
exponents z. More precisely, we only find continuous transitions for z < 1.46 or z > 6.98. 
The bounds are attained in the case of vanishingly small mass. 



By integrating the Gauss law ( 13 ) we obtain the charge density as 

Q = f * ^Lddr + i lim / * [Z($)F] = Q mcs . + Q frac . . (47) 

Here t\ and T2 are the radii between which a nonzero density of the charged fluid is present. 
These are the radii for which rh < £i\oc.( r )- For the mesonic phase r<i —> oo, as the charged 
fluid extends all the way into the IR. For fully fractionalized phases, in contrast, there are 
no radii for which the fluid is present. In general, given a numerical solution, we can read 
off the ratio Qfrac. /Q of fractionalized to total charge density. This ratio is plotted in figure 
[H] below as a function of the relevant UV coupling <pQ for the same values of {/3,m} as we 
used in the previous free energy plots. 

The plots exhibit two transitions, between mesonic and fractionalized and between par- 
tial and full fractionalization. As is apparent in the plots, the latter phase transition is 
softer than the former. The softer transition occurs as the width Ar = r<i — r\ occupied 
by the fermion fluid vanishes. The transition is very similar to that which occurs upon 
heating up the electron star without the dilatonic deformation [28, 29j. That transition 
was found to be third order. We refer the reader to those papers for more details, as this 
transition is not our main concern here. The third order nature of the transition can be 
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Figure 6: Ratio of fractionalized charge to total charge Q{ ra ,c./Q as a function of the relevant 
coupling for {rhj} = {0.1,20} (left) and {m,/3} = {0.5,10} (right). The dashed 
blue lines indicate a mesonic phase while the solid red lines indicate partially and fully 
fractionalized phases. The black dot denotes the location of the third order transition 
between partial and full fractionalization. The dashed line in the right hand plot indicates 
the jump in the ratio at the first order transition. 



heuristically understood as follows |28j . The difference in free energy due to the presence 
of the fluid is just the pressure of the fluid integrated over the radii at which it is present: 
ACl ~ Ar • p ~ Ar • 5fi^ . Here we defined 5(ii oc . = fi\ oc , — rh and used the expression (JtJ) 
for the pressure of the fermion fluid when the local chemical potential is only slightly bigger 
than the mass, as pertains close to the critical point where the fluid is about to disappear. 
Close to the transition one finds that Ar ~ ^P^ic ~ I A ~~ Acl 1 ^ 2 - Putting together the above 
formulae gives the third order transition Ail ~ \/l — fic\ 3 - 

By varying {/3, m}, it may be possible to achieve a first order transition directly from a 
mesonic to a fully fractionalized phase. This seems most likely for large masses. For masses 
bigger than the critical mass of figure [2] there is no Lifshitz solution to mediate between the 
mesonic and fractionalized phases. The transition is then necessarily first order. 

While Ofrac./O < 1 is a well defined bulk measure of the presence of mesino charge 
carriers in the bulk, it does not provide a useful field theory characterization of this fact. 
Two field theory descriptions of the same information can be obtained either by integrating 
the occupation numbers of the mesino fermion operator bilinear [481 [2"B] 

» = /^«^,*>, <«) 
or by summing up the volume of all the 2+1 dimensional Fermi surfaces in the boundary 
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mesino correlators [H [6] 



Q - Qfrac. = 2^ (27r) 2 Vn ' 



(49) 



The bulk WKB limit, that we have taken in order to use a fluid description for the fermions, 



implies that there are many Fermi surfaces in the sum (49). We can think of these as 



different excited mesino states. The WKB limit also implies that the tails of the fermion 
wavefunctions go to zero very quickly outside the boundary of the fluid, and so can be 
neglected for many purposes. Nonetheless, it is from these tails that we would be able to 



extract the expectation value in ( 48 ) 



Finally, let us return to the continuous fractionalization transition. The behavior of the 
ratio of fractionalized charge in figure [6] is suggestive of a second order phase transition, 
despite not being a symmetry breaking order parameter. To verify that this is the case, 
we can ask whether the charge susceptibility d 2 (l/djl 2 is discontinuous across the transi- 
tion. This is equivalent to the gradient of the charge density Q ~ d^l/djl jumping across 
the transition. The charge density can be extracted quite robustly from our numerics, as 
can, using the critical Lifshitz solution, the critical values of the ratios Q/ft 2 and 4>q/(m. 
Nonetheless, the question turns out to be delicate. In figure [7] below we plot data points 
very close to the phase transition, together with fits of the form y = ax + bx 2 . The slope a 
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Figure 7: Jump in the gradient of the charge density across the continuous fractionalization 
transition. Data points together with fit. 



in the fit jumps from ~ 0.269 on one side of the transition to ~ 0.253 on the other. This is 
a change of about 6%. We have tried to ensure that our numerics are sufficiently accurate 
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to pick out this small jump and hence conclude that the transition is indeed second order. 
We hope in the future it will be possible to rule out a higher order transition analytically. 

5 Final remarks 

The primary result of this paper is the construction of a simple holographic framework 
where the fraction of the charge density associated to gauge invariant Fermi surfaces could 
be tuned as a function of a relevant coupling in a quantum field theory. We exhibited 
phase transitions between 'mesonic' phases, wherein the Fermi surfaces account for all of 
the charge, to partially and fully fractionalized phases where they do not. At least for our 
models, this shows that neither the presence nor absence of gauge invariant Fermi surfaces 
is generic in holographic theories, it depends on relevant couplings in the theory. 

The phases we found broadly correspond with those anticipated in [T8] as FL, FL* and 
NFL phases. A feature of the particular model that we focused on is that it distinguished 
deconfmement of the glue sector from fractionalization of the charged sector. Our theories 
were gapless in all three phases, fractionalized or not. The class of models we have considered 
also admit gapped phases, and it is certainly of interest to study the interplay of confinement 
and Higgsing with fractionalization, in order to exhibit phase transitions into truly Fermi 
liquid phases [IHl [6] . 

For convenience we have worked in a coarse-grained fluid limit. This limit has the 
feature that there are many bulk Fermi surfaces. It has recently been demonstrated that 
it is technically feasible to go beyond this limit [6], and it will be important to generalize 
our discussion in that direction. That said, given that we start with a fractionalized, i.e. 
gauge-theoretic, description of the theory at large N, it is perhaps not surprising that once 
a single bound state forms, it is easy to form many further bound states. Whatever is the 
correct theory of quantum gravity for the universe in which we actually live, it must be 
compatible with the fact that we are surrounded by fluids. The fluid limit is natural in the 
bulk description and may yet have a deeper role to play in holography. This observation 
was also one of the motivations behind the earlier paper [27] . 

There is likely to be a bosonic analogue of the phase transitions we have discussed. In 
that case the mesonic phase will be superconducting rather than composed of Fermi surfaces 
|49j , and the onset of complete fractionalization will be accompanied by the recovery of the 
global U(l) symmetry. In fact, interesting phase transitions with this flavor have already 
been obtained in an Einstein-Maxwell-dilaton-charged scalar theory in [50] . 
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Throughout our phase diagram, the geometries become singular in the far interior, as 
we have discussed above. All of these singularities are well behaved in the sense that they 
can be regulated by heating up the system by a very small amount. Nonetheless, the 
singularities indicate that the true far IR of the system at zero temperature may revel new 
physics. In particular, it is possible that in the fractionalized phases e.g. higher derivative 
corrections will stabilize the dilaton at a highly curved AdS2 x M 2 IR fixed point, along the 
lines of [51], or alternatively gap the theory. 

The fractionalization transitions that we have found do not realize one of the key features 
desired in |22j . Namely, the fractionalization proceeds by the Fermi surfaces successively 
vanishing as their Fermi momenta kp — > 0, as the local bulk chemical potential drops below 
the fermion mass, in contrast to the residues vanishing, Z — > 0, with the Fermi surface size 
remaining finite. While it is possible to obtain holographic Fermi surfaces with vanishing 
residue given an AdS2 x M 2 near horizon geometry [TU] , those Fermi surfaces are not in fact 
themselves fractionalized (although they do interact with a fractionalized sector described 
by the extremal horizon), precisely because they appear as singularities in gauge invariant 
mesino Green's functions. To obtain a holographic fractionalization transition in the spirit 
of [22] one must take the theory into a Higgsed phase. 

Finally, it is of interest to define an order parameter, a generalization perhaps of the 
Polyakov loop, that is capable of distinguishing fractionalized and mesonic phases. Our 
phase transitions give a test bed in which to construct and understand such an order param- 
eter holographically. Once defined as a field theoretical quantity, such an order parameter 
should be useful beyond holography. 
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